Density Matrix Renormalization Group Study 
of the Haldane Phase in Random One-Dimensional Antiferromagnets 
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It is conjectured that the Haldane phase of the S = 
1 antiferromagnetic Heisenberg chain and the S = 1/2 
ferromagnetic-antiferromagnetic alternating Heisenberg chain 
is stable against any strength of randomness, because of im- 
posed breakdown of translational symmetry. This conjecture 
is confirmed by the density matrix renormalization group cal- 
culation of the string order parameter and the energy gap 
distribution. 
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In the recent studies of quantum many body problem, 
the ground state properties of the random quantum spin 
systems have been attracting a renewed interest |H-E5|. 
Among them, the effect of randomness on the spin gap 
state of quantum spin chains has been extensively studied 
theoretically and experimentally 

The real space renormalization group (RSRG) method 
has been often used for the study of random quantum 
spin chains. Using this method, it has been exactly 
proved that the ground state of the 5 — 1/2 random 
antiferromagnetic Heisenberg chain (RAHC) is the ran- 
dom singlet(RS) state |l|-§] irrespective of the strength of 
randomness. Hyman et al. ]5| have applied this method 
to the S — 1/2 dimerized RAHC and have shown that 
the dimerization is relevant to the RS phase. They con- 
cluded that the ground state of this model is the random 
dimer (RD) phase in which the string long range order 
survives even in the presence of randomness ]5|,§| . These 
results are numerically confirmed using the density ma- 
trix renormalization group (DMRG) method . 

The effect of randomness on the Haldane phase is also 
studied by Boechat and coworkers (^|J^] and Hyman and 
Yang H using the RSRG method for the original model 
and the low energy effective model, respectively. These 
authors predicted the possibility of the RS phase for 
strong enough randomness. This problem has been fur- 
ther studied by Monthus and coworkers using the numer- 
ical analysis of the RSRG equation for the square distri- 
bution of exchange coupling Jl(|. They predicted that 
the Haldane-RS phase transition takes place at a finite 
critical strength of randomness. In the finite neighbour- 
hood of the critical point, the Haldane phase belongs to 
the Griffith phase with finite dynamical exponent z > 1. 
Hereafter this phase is called the random Haldane (RH) 
phase. On the other hand, Nishiyama Jll]] has carried out 



the exact diagonalization study of the 5 = 1 RAHC. He 
observed that the Haldane phase is quite robust against 
randomness and the string order remains finite unless the 
bond strength is distributed down to zero. He also car- 
ried out the quantum Monte Carlo simulation Jl2] ] and 
found no random singlet phase even for strong random- 
ness. On the contrary, the quantum Monte Carlo simu- 
lation by Todo et al. jl3| suggested the presence of the 
RS phase for strong enough randomness. 

In the absence of randomness, the present author has 
given a physical picture of the Haldane phase as the lim- 
iting case of the 5 = 1/2 Heisenberg chain with bond 
alternation in which the exchange coupling takes two dif- 
ferent values J and Jp alternatingly jn| . In the extreme 
case of Jp — > — oo, this system tends to the 5 = 1 antifer- 
romagnetic Heisenberg chain. The string order remains 
long ranged over the whole range — oo < Jp < J and only 
vanishes at J = Jp. The perfect string order is realized 
for Jp = 0. As discussed by Hyman et al this is the 
direct consequence of the imposed breakdown of transla- 
tional symmetry. Because the randomness cannot recover 
the translational symmetry, the string order is expected 
to remain finite over the whole range — oo < Jp < J 
for any strength of randomness. Therefore we may safely 
conjecture that the Haldane phase of the 5 = 1 RAHC 
should also remain stable for any strength of randomness. 

In the following, we confirm this conjecture using the 
DMRG method 0,[l7j which allows the calculation of the 
ground state and low energy properties of large systems 
with high accuracy. We use the algorithm introduced 
in ref. Q. This method has been successfully applied 
to the spin- 1/2 RAHC and weakly dimerized spin-1/2 
RAHC in which the system is gapless or has very small 
gap in the absence of randomness. Namely, in these sys- 
tems the characteristic energy scale of the regular sys- 
tem is much smaller than the strength of randomness 
even for weak randomness. Compared to these examples, 
present model is less dangerous because the regular sys- 
tem has a finite gap and the characteristic energy scale 
of the regular system is comparable to the strength of 
randomness even in the worst case. We investigate not 
only the 5 = 1 RAHC but also the 5 = 1/2 random 
ferromagnetic-antiferromagnetic alternating Heisenberg 
chain (RFAHC) which interpolates the 5 = 1/2 dimer- 
ized RAHC and the 5 = 1 RAHC. 

The Hamiltionan of the 5 = 1/2 RFAHC is given by, 
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ti = ^ 2JjS , 2 i-iS , 2i + 2J F S , 2iS , 2 i + i, | Si |= 1/2, (1) 



i=l 



where Jp = const, and J^'s are distributed randomly 
with probability distribution, 



P(Ji) 



l/W for 1 - W/2 < Ji < 1 + W/2 
otherwise. 



(2) 



The width W of the distribution represents the strength 
of randomness. The maximum randomness is defined by 
W = 2, because the ferromagnetic bonds appear among 
Ji's for W > 2. It should be noted that the appearance 
of the random ferromagnetic bonds can drive the system 
to the completely different fixed point called large spin 
phase JTsfl . Although the crossover between the random 
Haldane phase and the large spin phase is an interest- 
ing issue, we leave this problem outside the scope of the 
present study. 

The ground state of the regular counterpart of this 
model (W — 0) is the Haldane phase with long range 
string order defined by O s i T = lim^jv-xx, O s ^ T (l; N) El- 
Here O s i r (l; N) is the string correlation function in the 
chain of length N defined only for odd I as, 



{2i+;+i i 
fc=2i+l J 



(3) 



where < ... > denotes the ground state average. In the 
presence of randomness, the string order is defined as 
the sample average of O s t r . In the limit Jp — * — oo, 
the string order parameter (^) reduces to the one for the 
5 = 1 chain 

First, we calculate s t r for the S = 1/2 RFAHC and 
5 = 1 RAHC using the DMRG method. The calcula- 
tion is performed with open boundary condition. For the 
5 = 1/2 RFAHC, the bonds at the both ends of the chain 
are chosen to be antiferromagnetic to avoid the quasi- 
degeneracy of the ground state. The two boundary spins 
are not counted in the number of spins 27V to keep the 
consistency with the 5 = 1 chain (see below). The av- 
erage is taken over 200 samples with N < 29 (58 spins). 
The string order for the finite system O s ^ T (N) is esti- 
mated from O s i T (l; N) averaged over 6 values of I around 
I = N/2. The maximum number m of the states kept in 
each step is 100. Similar calculation is also performed 
for the 5=1 RAHC. In this case, the additional spins 
with 5=1/2 are added at the both ends of the chain 
to remove the quasi-degeneracy as proposed by White 
and Huse for the regular chain J2lJ . The average is taken 
over 400 samples with N < 42 where N is the number of 
5=1 spins. In this case, we take m = 80 and O s ^ T (N) 
is estimated from O s ^ T (l;N) averaged over 12 values of 
I around I = N/2. We have confirmed that these val- 
ues of m are large enough from the m-dependence of the 
obtained values of O s ^ T (N). 




FIG. 1. The size dependence of the string order parame- 
ter O s ^ T (N) plotted against 1/N. The solid lines are fit by 
the formula (6). The points with different values of W are 
depicted by the same symbols for W = 1.6, 1.8, 1.9 and 2.0 
from top to down. 



0.5 







□ 


O m 


□ 


W=0 » regular 


• U 

□ 




1.6 • 






1.8 o 

1.9 ■ 






2.0 = 







0.1 



0.2 



1/Mf 



FIG. 2. The Jp-dependence of the string order parameter 
O str for W = 1.6, 1.8, 1.9 and 2.0 plotted against 1/ | J F |. 

Figure [j] shows the size dependence of the string order. 
Typical sizes of the error bars estimated from the statisti- 
cal flucuation among samples are less than the size of the 
symbols unless they are explicitly shown in the figures. 
The extrapolation is made under the assumption 



OstrW =s Ostr + CN^ exp(-A/£). 



(4) 



where C and £ are the constants to be determined by 
fitting. The exponent 2r\ characterize the size dependence 
of the string order parameter at the RH-RS critical point 
where £ should diverge as O s ^ T (N) ~ N~ 2ri . This value is 
estimated as follows: According to Monthus et al. Jl(|, at 
the critical point, O s ^ r (A f ) behaves as r~ 2 ( 3 ~*) with <j) — 
y/E while the logarithmic energy scale T varies with the 
system size as V ~ A 1 / 3 . Therefore O s ^ r (7V) should scale 
as Ar- 2 (3-0)/3 ^ ^y-0.5092 at tlie cr i t } ca j pomt resulting 

in 2i] = 0.5092. It should be noted that the low energy 
effective model of Hyman and Yang || also applies for the 
5 = 1/2 RFAHC with finite Jp < by construction. The 
extrapolated values of O s tr are plotted against 1/ | Jp | 
in Fig. |[ The string order is perfect at Jp = 0, where 
the ground state is a simple assembly of local singlets 
and should decrease with the increase of | Jp | . This 
behavior is clearly seen in Fig |[ In this extrapolation 
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scheme, the string long range order remains finite even 
at W = 2 for both 5 = 1/2 RFAHC and 5 = 1 RAHC. 




FIG. 3. The size dependence of the string order parameter 
O s tr(N) for the 5" = 1 RAHC plotted against jV" - 5092 . 




FIG. 4. The IV-dependence of the string order parameter 
O str for the S = 1 RAHC. 

To check our extrapolation scheme, we also made the 
extrapolation assuming the size dependence N~ 2ri ex- 
pected at the RH-RS critical point in Fig. || for 5 = 1 
RAHC. In the RH phase, the extrapolated values thus 
obtained can be understood as the lower bound. For 
W > 1.9, it is clear that the extrapolated values remain 
definitely positive. Therefore the extrapolation using Eq. 
(0) is appropriate in this region rather than the power law 
extrapolation. For W = 2, the extrapolated value is very 
small but still positive (~ 0.018). 

Even if we do not rely on the values for W = 2 ex- 
trapolated using Eq. (0), we can convince ourselves the 
stability of the RH phase at W = 2 by the following ar- 
gument. In Fig. ^, we plot O s t r extrapolated using Eq. 
@ for 5 = 1 RAHC against 2 - W for W < 1.9. If we 
assume the critical behavior O s ^ r ~ (W c — W) 1 - 173 pre- 
dicted by Monthus et al. [^0) , it is highly unlikely that the 
string order disappears at finite values of W c less than 2. 

Especially, this excludes the possibility W c ~ 1.485 
predicted by Monthus et al. jl0|. In general, it is not 
surprising that the RSRG method gives incorrect value 
for the critical point even if it gives correct values for 
the critical exponents, because the RSRG transforma- 
tion is not exact at the initial stage of renormalization. 
Furthermore, Monthus et al. Jl(| have neglected the ef- 
fective ferromagnetic coupling between the next nearest 
neighbour interaction which appear after decimation of 
two 5 = 1 spins. (See the discussion following eq. (2.20) 
of ref. Jk}.) The neglect of this term is equivalent to the 



introduction of the antiferromagnetic next nearest neigh- 
bour interaction as a counter term in the bare interac- 
tion. In terms of the RFAHC, such interaction makes the 
distinction between the even and odd bonds meaningless 
and can recover the translational symmetry leading to 
the destruction of the string order erroneously. 




FIG. 5. The iV-dependence of < ln(l/A) > (open symbols) 
and a (filled symbols) of the S = 1 RAHC plotted against 
InN for W = 1.8, 1.9 and 2.0. The horizontal arrows indicate 
the values of 1/Po estimated from d < ln(l/A) > /dlnN for 
W — 2, 1.9 and 1.8 from top to bottom. 

To further confirm the stability of the RH phase at the 
most dangerous point Jp — > — oo and W<2, we calculate 
the energy gap distribution for the 5 = 1 RAHC. In the 
RH phase, the fixed point distribution of the energy gap 
A is given by P(x) = Po exp(— Pqx) where x = In(f2/A), 
fl is the energy cut-off and Pq is the nonuniversal con- 
stant ||. For the finite size systems, ti scales as N^ 1 ^ 
||. Therefore the dynamical exponent z is given by 
z = 1/Pq. Furthermore, this distribution implies 



< ln(l/A) >= 



Pn" 1 hi N 



const. 



(5) 



o = y/< (In A- < In A >) 2 > = P^ 1 = z (6) 

for N >> 1. In Fig. [|, we pl ot < hi(l/A) > and a 
against IniV. The error bars are estimated from the sta- 
tistical flucutation among samples. The average is taken 
over more than 100 samples with N < 50. For the most 
random case W = 2, the average is taken over 219 sam- 
ples. In this case, we have taken m = 100 in most cases. 
For the confirmation of the accuracy, however, we recal- 
culated with m = 150 for the samples with very small 
gap (less than 10~ 3 ) but the difference was negligible. 
Actually, the latter data are also plotted in Fig. ^ for 
W = 2. But they are almost covered by the data with 
m = 100 and are invisible in Fig. ^|. Therefore we may 
safely neglect the m-dependence for less dangerous case 
W < 2. 

It is evident that < ln(l/A) > behaves linearly with 
InJV and a tends to a constant value as N — > oo. We 
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can estimate the values of z(= l/-Po) from the gradient 
d < ln(l/A) > /dlnN for large N. These values are 
indicated by the arrows in Fig. || and they are consistent 
with those estimated from a for large N within the error 
bars. On the other hand, Figs. ra(a) and (b) show the 
plot of < ln(l/A) > against A 1 ™ and A 1 / 2 which are 
the expected size dependence at the RH-RS critical point 
and within the RS phase, respectively |T(|. Both plots 
are less linear compared to Fig. ||. These results confirm 
that the ground state of the 5=1 RAHC remains in the 
RH phase down to W = 2. 



order and the energy gap distribution. 

The numerical calculations have been performed using 
the FACOM VPP500 at the Supercomputer Center, In- 
stitute for Solid State Physics, University of Tokyo. The 
author thanks H. Takayama and S. Todo for useful dis- 
cussion and comments. This work is partially supported 
by the Grant-in- Aid for Scientific Research from the Min- 
istry of Education, Science, Sports and Culture. 
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FIG. 6. The TV-dependence of < ln(l/A) > of the S = 1 
RAHC plotted (a) N 1/3 and (b) N 1/2 for W = 1.8, 1.9 and 
2.0. 

It should be also noted that the finite size effect be- 
comes serious only if one hopes to conclude the presence 
of the RS phase. Even in the RH phase, the string or- 
der or the gap distribution might behave RS-like if the 
system size is not enough. Actually, the authors of ref. 
||l0f needed extremely large number of spins to conclude 
that their calculation leads to the RS phase. But the RS 
phase can never behave RH-like by the finite size effect 
because the RS phase has divergent correlation length. 
Therefore, it is relatively easy to exclude the possibility 
of RS phase if the deviation from the RS-like behavior 
is already observed for relatively small systems, which is 
the case of the present calculation. 

In summary, it is conjectured that the Haldane phase 
of the 5 = 1 RAHC and the 5=1/2 RFAHC is stable 
against any strength of randomness, because of the im- 
posed breakdown of translational symmetry. This conjec- 
ture is confirmed by the DMRG calculation of the string 
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